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Abstract 

PHENIX reports on electromagnetic and hadronic observables in large data sets of p+p, d+Au 
and Au+Au collisions at various cms energies. Initial state effects in cold nuclear matter are 
quantified by centrality dependent n°, tj, reconstructed jets and if// measurements. Using the 
first run of the new EBIS ion source at RHIC, we report first results for particle flow (vi and 
v'2) and quarkonium production in U+U and Cu+Au collisions. Hot matter created in Au+Au is 
characterized using event-plane dependent HBT and dielectrons. Parton-medium interactions are 
investigated using high pj single hadrons, y-hadron correlations and heavy flavor decay electrons 
identified with the newly installed VTX detector. 



1. Introduction 

RHIC experiments have demonstrated that quark-gluon plasma (QGP) is created in Au+Au 
collisions. Subsequent experimental goals focus on characterizing the QGP, and exploring the 
phase transition and a possible critical point by varying the colliding nuclei and collision energy. 
We present the latest results from the PHENIX experiment at RHIC. The results are categorized 
according to time within the collision: initial state effects, hot matter dynamics, and parton- 
medium interactions. 



2. The baseline for QGP properties 

Initial state effects on the primordial hard scatterings which produce hard probes of the hot 
dense medium reflect parton distribution functions in nucleons and nuclei, as well as possible 
energy loss or bound state dissociation in cold matter, p+p collisions provide the baseline, while 
d+Au collisions reflect additional cold nuclear matter (CNM) effects arising from the presence of 
a nucleus in the collision. Direct photon spectra in d+Au collisions show little or no modification 
compared to expectations from p+p collisions, though an isospin effect is consistent with the 
data IHI2]. On the other hand, hadronic probes display centrality-dependent CNM effects. 

In Fig.Q] the nuclear modification factors (RdA) for n°, r\ and fully-reconstructed jets in d+Au 
collisions at y/s?w=200 GeV are shown [3]. These three probes result from different analysis 
techniques, but show remarkable consistency. While a slight suppression is found in the most 
central collisions, there is significant enhancement at high pj in the peripheral collisions; such 
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Figure 1: Nuclear modification factors (Rja) for w \ *7 and fully-reconstructed jets in d+Au collisions. 



a feature is not predicted by available models. We note that centrality is determined in the data 
by the particle multiplicity in the beam-beam counter (BBC) at 3.1 < \rj\ <3.9, while the models 
use impact parameters. Therefore, the centrality classes between data and models may differ 
somewhat, particularly for events with high pj hadrons or jets (pj >10GeV). 

Quarkonia offer key probes of initial state effects. The high statistics 2008 data enable mea- 
surement of \pi production as a function of collision centrality in d+Au at ys/wv=200 GeV. The 
dielectron mass spectrum in the J/if/ and ipt region is shown in Fig.|2](|3t]. The ratio of i////(J/t//) is 



d+Au 0-88% Centrality 




M„.„. [GeW] 



Figure 2: iftr peak seen in dielectron mass spectra in 
minimum bias d+Au collisions. 
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Figure 3: Rja of ipr and J /if/ in d+Au collisions as a 
function of centrality. ipi is more suppressed than J/ift. 



~ 0.8 %, smaller than in p+p (2.1 %). The centrality dependence of the ratio is shown in Fig. [3] 
The yield of J/ijj is somewhat suppressed in d+Au collisions compared to p+p, but t/r/ is more 
heavily suppressed. The relative suppression at FNAL fixed target energy can be reasonably well 
explained by considering the time the precursor cc state spends in the nucleus. However, this 
explanation completely fails at RHIC energy. 

3. New collision systems for QGP characterization 

In 2012, first data were taken with the EBIS ion source, providing new ion species at RHIC. 
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3.1. U+U collisions 

Uranium is a highly deformed nucleus. Therefore, selecting collisions where two uranium 
nuclei are aligning with their long axes head-on (tip-tip) should access higher energy density at 
RHIC energy. PHENIX collected data in U+U collisions at tJsnn— 193 GeV. First results show 
that the Bjorken energy density for 0-2 % centrality (tip-tip enriched events) is ~30 % higher than 
the most central Au+Au collisions at 200GeV 04J]. The elliptic flow of identified charged pions 
and protons in U+U collisions are compared to that in Au+Au at 0-10% centrality in Fig.|4]|5]. 
For pions in U+U, V2 is flatter with pj for pj <1 GeV, reminiscent of observations in Pb+Pb at 
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Figure 4: Elliptic flow of charged pions and protons in (a) 
Au+Au and (b) U+U collisions at 0-10 % centrality. 



Figure 5: Elliptic flow of charged pions and protons in 
(a) 0-2 % and (b) 6-10 % centrality in U+U collisions. 



LHC (a]. This can be interpreted as a consequence of very strong radial flow for higher energy 
density matter. The finer centrality classes shown in Fig. [5] show that this flattening is most 
pronounced in 0-2 %, where tip-tip events represent a larger fraction of the total, and is nearly 
nonexistent for 6-10 % centrality. It is intriguing to see this phenomenon at RHIC, as the energy 
density reached at LHC is ~3 times higher than at RHIC. 

3.2. Cu+Au collisions 
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Figure 6: i'i of charged pions in Cu+Au collisions. 



Figure 7: V2 of charged pions and protons in Cu+Au 
collisions. 



In Cu+Au collisions, one can select events in which the Cu is fully contained within the Au 
nucleus. Furthermore, the pressure gradient developed should be asymmetric, resulting in an 
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intrinsic triangularity in the collective motion. This triangularity generates a V3 that is geometric, 
rather than fluctuation driven. PHENIX has already analyzed ~20% of the 2012 Cu+Au data. 
The shower-max detector in the zero-degree calorimeters is used to determine the event plane 
(reaction plane *fi), allowing a measurement of vi. We defined the Au-spectator going side as 
positive 

Fig.|6]and|7]show vi and vi for identified pions and protons. A sizable positive \\ is observed; 
the sign and magnitude of v\ differ from predictions of the the AMPT model (HIJING+parton 
cascade) |7|], although AMPT describes symmetric collision systems very well. 

Utilizing the muon arms at 1.2< \r\\ <2.2, PHENIX measured J/t// spectra in Cu+Au colli- 
sions via the di-muon decay channel. This allowed determination of J/t/i Raa in Cu+Au, and 
comparison to the the J/t// suppression observed in Cu+Cu and Au+Au collisions. The results 
and comparison are shown in Fig. [8] [8] . The suppression of the J/t// yield as a function of N part 
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Figure 8: Raa of J/ iff in Cu+Au, Cu+Cu and Au+Au collisions. 

on the Au-going side fits well with the N part dependence in Au+Au and Cu+Cu collisions. How- 
ever, on the Cu-going side, the suppression is significantly stronger. The difference of the nuclear 
PDFs in Cu and Au nuclei partly explains this trend, but additional final state effects may need 
to be considered in order to understand the magnitude on the Cu-going side. 

4. Characterizing hot dense matter by new probes 

Two boson correlations as pioneered by Hanbury-Brown and Twiss (colloquially known as 
HBT correlations in our field) provide information on the space-time evolution of the particle 
emission source. PHENIX has measured the angular dependence of HBT radii with respect to 
the (second order) event plane {R 2 (A(f>) = R 2 Q + 2 x I.R 2 n cos(nA(f>)). This analysis shows that the 
eccentricity, defined as €hbt = obtained from the Kaon HBT correlations is consistent 

with the initial eccentricities. This implies early freezeout of kaons from the expanding hadronic 
phase late in the collision. 

PHENIX studied the HBT radii with respect to the 3rd order event plane for the first time. 
Fig. [5] shows R s and R for charged pions as a function of the emission angle with respect to the 
2nd and 3rd order event plane in 0-10 % in Au+Au collisions [9]. The average radii are 10 and 5 
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Figure 9: Angular dependence of HBT radii (R s and R ) with respect to 2nd (¥2) and 3rd (¥3) order event plane in 
0-10% Au+Au collsion at 7^=200 GeV. 



for ¥2 and ¥3, respectively. While R s shows only a weak dependence on orientation with respect 
to the event plane, R indicates a clear temporal variation as a function of angle for both the 
¥2 and ¥3 event planes. Moreover, the magnitude of the oscillation for ¥3 is almost the same 
as ¥2. Since R a includes the emission duration in addition to the spatial homogeneity length 
in the transverse direction, this measurement reflects the space-time evolution of ellipticity and 
triangularity of the collision system. 

Turning to thermal radiation, we report on a unique and penetrating probe of thermodynamic 
quantities. Measuring thermal radiation gives access, for example, to the temperature of the hot 
dense medium produced in the collision. Using the Hadron Blind Detector (HBD) installed in 
2010, PHENIX performed a new measurement of thermal di-electron radiation. In all bins from 
peripheral to semi-central, the new result is consistent with that from the 2004 data. The most 
central bin is still being analyzed. For further details, see IflQII . 



5. Detailed study of energy loss of partons 

5.1. y-h correlation 



5<pi f <9GeV 

[ 1 

H iH 


7c x 0.5 < < 7 GeV/c i 

. \A$ > ir/2 
*, . \A<$,\ > 2ji/3 -_ 

Ti t ffll > 5lt/6 ~~ 




_t ^ 0-40% Au+Au = 

ffl PH#ENIX i 

preliminary = 



0.5 1 1.5 2 2.5 3 „ 3.5 



Figure 10: Ratio of per-trigger away-side hadron yield in 200 GeV Au+Au and p+p y-hadron measurements as a function 
off. The legend identifies the three different integrating angular ranges. 



y-hadron correlations are often considered a "golden channel" to evaluate the energy loss 
of partons in the medium. Since direct photons reflect the momentum of the original scattered 
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partons, they tag the initial momentum of the scattered parton that traverses the medium. Conse- 
quently one can deduce the energy loss via bpj = pj(y) - pr(hadrons). PHENIX has measured 
the hadron yield opposing direct photons as a function of £(= -ln(zr)) for the 2007 and 2010 
Au+Au data sets together. The results are shown in Fig. [TOl 111 ill . As the angular range in which 
the away-side hadrons are integrated is increased from \6<f>\ > 5n/6 to \6<f>\ > n/2, the per-trigger 
away-side yield increases compared to p+p collisions. The increase is in the larger £ or small zt 
region, implying that the softer away-side jet fragments are distributed over a wider angle than 
fragments carrying a large fraction of the jet energy. This is expected from medium-enhanced 
splitting, and is qualitatively consistent with the angular broadening observed in hadron-hadron 
correlations at RHIC and the jet shape broadening observed at LHC. 

5.2. High pt n° to determine fractional momentum shift 

Using the large 2007 Au+Au data set at yjs^N = 200 GeV and a new algorithm to correct for 
shower merging, PHENIX extended the pj range of 7r° spectra 111 211 . Fig. Q~T] compares n° Raa in 
200 GeV Au+Au collisions from RHIC and charged hadron Raa measured in 2.76 TeV Pb+Pb 
at the LHC (ALICE experiment) lfl3ll . For both centralities and over the entire pj range shown, 
the two data sets are rather similar. However, this does not necessarily mean that the energy 
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Figure 11: Comparison of Raa of tP from PHENIX at RHIC and charged hadrons from ALICE experiments at LHC. 
Two data sets from ALICE experiments correspond to two different p+p references. 



loss of the partons is the same, as the spectral shapes differ. PHENIX determines the average 
fractional momentum shift (S i oss ) of high pj n to estimate the average fractional energy loss 
of the initial parton. S i oss is defined as dpr/pr, where 6pj is the difference of the momentum 
in p+p collisions (f>T,pp) and Au+Au (/?t,AuAu); pr in the denominator is pj tVP . The calculation 
method is schematically depicted in Fig. Q~2] First, the jt° cross-section in p+p (f(pj)) is scaled 
by Taa corresponding to the centrality selection of the Au+Au data (g(pj)). Next, the scaled 
p+p cross-section (7\a./Xpt)) is fit with a power-law function (h(pj)). Third, a scaled p+p point 
closest in yield (p' T ) to the Au+Au point of interest (/?t,AuAu) is shifted using the fit function as 

WCPT.pp) = h(p T , pp )/h(p' T pp ) x WCPx.pp), (1) 
where /?T,pp is chosen to satisfy the relation 



7aa/(PT,pp) = g(/?T,AuAu)- 
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(2) 



(2) Move scaled p+p 
data, along fit 



(3) Calculate 5p T 
□ = pT(p+p) - pT(Au+Au) 



p+p data ' 
and fit tunc. *' 



10 11 12 13 14 15 

p T (GeV/c) 



-> H 0.45 

Q. 

a K 0.4 

CO 

X 0.35 
«~ 

0.3 
0.25 

0.2 
0.15 

0.1 
0.05 



Pb+Pb 0-5%, S(global)=0.3% 
Au+Au 0-5%, 5(global)=1 .0% 
Pb+Pb 70-80%, fi(global)=0.7% 
Au+Au 70-80%, 8(global)=2.9% 



B 



16 18 20 
p T (p+p) (GeV/c) 



Figure 12: Demonstration of calculating the frac- 
tional momentum loss (S/ oss ). 



Figure 13: 5 /„.„ for PHENIX n° and ALICE charged 
hadrons. 



The 6pr is calculated as /?T,pp - /?t,AuAu- For obtaining S i oss , the 6pj is divided by the pj,pp- 
We show the comparison of S i oss for PHENIX 7r°'s in Au+Au and ALICE charged hadrons in 
Pb+Pb in Fig. [T3] We took the corresponding Pb+Pb and p+p data for black points in Fig. QT| 
for calculating the S i oss for ALICE charged hadrons. The S i oss for the ALICE charged hadrons 
is found to be ~25 % higher than that for the PHENIX 7r°'s. PHENIX has also measured S i oss in 
Au+Au at 39 and 62.4 GeV cms energy. The change in S i oss reaches approximately a factor of 4 
from 200 GeV to 39 GeV H, [3 . 



6. Heavy flavor electrons 

The suppression of non-photonic electron yields was found to be nearly as large as that of 
7r°'s. The question has been whether the suppression arises from the suppression of charm and/or 
bottom quarks. Theoretical predictions about this differ greatly. Consequently, PHENIX con- 
structed a micro vertex detector to separate electrons from charm and bottom decays. This detec- 
tor, the VTX, consists of 2 layers each of strip-pixel and pixel silicon detectors. The measurement 
of the distance of closest approach (DCA) of electrons to the collision vertex allows separation 
according to the parent hadrons. Fig. [14]shows a DCA distribution measured in the VTX detector 
in p+p collisions. The distribution was deconvoluted to determine the relative contributions of 
the heavy flavor parent hadrons. Fig. [15] shows the ratio of electron yields from b-quarks, and 
b + c quarks for p+p collisions QSQ . The points with large error bars are from Kn correlation 
measurements 111 611 . Our new result is consistent with both the Kn correlation analysis and with 
FONLL calculations. We also analyzed electrons in Au+Au collisions to obtain the Raa for elec- 
trons from charm and bottom quarks [8]. The analysis assumed that the parent hadrons (e.g. D, 
B) follow the pj distributions given by PYTHIA. The results demonstrate that the Au+Au data 
are inconsistent with these input assumptions unless there is also a large suppression of electrons 
from bottom across the measured p e T range. However, this large suppression implies a change in 
the parent hadron pj distributions, which results in changes in the electron DCA distributions. 
We are actively working on evaluation of the uncertainties caused by this fact. 
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Figure 14: DCA distribution in p+p collisions for 2< 
p e T <2.5 GeV/c, and its deconvolution to contributions 
from parent hadrons. 



Figure 15: Ratio of electron yields from bottom 
quarks to bottom and charm quarks in p+p colli- 
sions. 



7. Summary 

PHENIX has reported new findings in collisions at various energies and species available at 
RHIC. There is a significant enhancement of high pj hadrons and jets in the peripheral d+Au 
collisions. The yield of if// is more heavily suppressed than that of J /iff in d+Au collisions. A 
strong radial flow is seen in tip-tip enriched 0-2% U+U collisions. A sizable positive vi is 
observed in Cu+Au collisions. In the same system, the J /if/ yield measured in 1.2< I77I <2.2 
as a function of N parr shows larger suppression in the Cu-going side than the Au-going side. 
The space-time evolution of the ellipticity and triangularity of the Au+Au collision system is 
observed in angular dependent HBT measurement. PHENIX has seen an angular broadening 
of the away-side soft particles in y-h measurement, similar to h-h correlations and jet shape 
measurements from RHIC and LHC. The high pj hadron spectra showed that the fractional 
momentum shift under the presence of hot dense matter is ~25 % higher at LHC compared to the 
RHIC top energy. Identification of electrons from bottom and charm quarks has been successful 
at RHIC using the VTX detector. The ratio of (b — > e)/(b, c — > e) for p+p collisions is consistent 
with the previous measurement and well described by a FONLL calculation. 
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